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fvq ■ Abstract. The Very High Energy (VHE) 7-ray emission of the BL Lacertae objects 

Markarian 501 and Markarian 421 has been observed by the CAT Imaging Atmospheric 

Cherenkov Telescope in 1997 and 1998. The spectrum extraction method is presented, 

^ I and the spectral properties of both sources are compared in different activity states. 

Theoretical implications for jet astrophysics are discussed. 
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SPECTRAL ANALYSIS METHOD 

The CAT (Cherenkov Array at Themis) telescope records Cherenkov flashes due 
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p i' to VHE atmospheric showers through its 17. 8m^ mirror. The fine grain of its cam- 

Q ! era, combined with fast electronics, allows a relatively low 7-ray detection threshold 

^ I energy of 250 GeV (at Zenith) as well as an accurate analysis of the resulting images. 

52 ' The experiment is fully described elsewhere [3]. The analysis method is based on 

the fit of each individual shower image with theoretical mean 7-ray images coming 
from a semi-analytical model [7] . This permits the extraction of a 7-ray signal from 
r> I the background due to cosmic-ray induced atmospheric showers; the optimization 

c^ ■ of the event selection using the Crab nebula data gives a hadron-rejection factor 

of ~200 for a 7-ray efficiency of ~40%. The fit also yields the energy of each 
shower in the 7 hypothesis, with a resolution of ~20%, independent of energy, 
when restricting to showers with a small impact parameter (< 130 m). However, 
the spectrum reconstruction procedure presented below involves the exact energy- 
resolution function "${£ ^ E, cosO) ^ in order to use the entire available statistics. 
These functions have been determined by detailed Monte-Carlo simulations of the 
telescope response, as has the effective detection area ^eff, which includes the effect 
of event-selection efficiency. The simulations have been checked and calibrated on 
the basis of several observables by using muons rings and the nearly-pure 7-ray 
signal from the highest fiare of Mrk 501 in 1997 [9]. 



Probability density, for fixed zenith angle 9 and real energy E, to get an estimated energy E. 



With typical statistics of ~1000 7-ray events and a signal-to-background ratio 
of 0.4 (as obtained on the Crab nebula), a spectrum can be determined with rea- 
sonable accuracy as follows. First we define a set {Aj^} = {[i?™"^, i?™^^]}j^=i^„^ 
of energy bins; since the threshold energy of the telescope increases with zenith 
angle, each of these is divided into a subset {Aj^^j^} = {[0f^™,9f^^^]}i^=i^n^{ie) of 
zenith angle bins ^ with a width of 6{cos6) = 0.05. Then, for all ON and OFF- 
source runs ^, the number of events passing the selection cuts is determined within 
each Ajg j^ 2D-bin. Finally, assuming a given spectral shape, a maximum-likelihood 
estimation of the spectral parameters is performed, which takes into account the 
effective detection area and the energy-resolution function. At this stage, two hy- 
potheses Ti^^ and T-f^^ are successively considered, which are a simple power law 

(f)P\E) = 0o^TeV and a curved shape 0^^(E) = 0o£'Tev^^'°'^'"^''°''\ respectively. 
Power laws often account for spectral properties in high energy astrophysics, at 
least because they often constitute a good approximation when considered over a 
restricted energy range. On the other hand, a curved shape is suggested by general 
considerations on emission processes within blazar jets (see the next sections); the 
latter parametrization, previously used by the Whipple group [11], corresponds to 
a parabolic law in a log(0) vs. log{E) representation and allows simple compar- 
isons. The relevance of TY^' with respect to Ti'^^ is estimated through the likelihood 
ratio of the two hypotheses, which is defined as A = — 2 x log(^): it behaves 
(asymptotically) like a x^ with one d.o.f. and permits a search for the presence 
of spectral curvature. All spectral curves shown in the next sections come from 
the best set of fitted parameters obtained by this procedure. It is also customary 
to superimpose on them "experimental" points {El"^^ ^(jf^^}i^=i^ne.; which are only 
indicative of the statistics used in each Aj^ energy bin. The definition of these 
points is the following: if we note (^"^^^{E) as the best fitted spectrum under the 
hypothesis 0^' or cjf'' according to A, N°^l^ as the number of 7-ray events observed 

in the A^^ i^ 2D-bin, and a^, i, as its error, and Nf'f = / J' dE / dE 0''""VE) x 

J _B™'" J 

Acs{E, {cos 9) i^)x'^{E -^ E, {cos9)i^) as the corresponding predicted number, then 
E^^'^ is the unique energy within the bin for which (p^^^^lE^^'^) = {(f)^^^^)^ holds ^, 

nz(ie) 
iz=l 

noted that this way of representing spectra cannot replace the likelihood method, 
which provides the only relevant physical results, i.e. the values of the most proba- 
ble spectral parameters as well as their errors and covariance matrix. In particular, 
one should not use the {E'^^^, 0^^^} points in any kind of minimization, since the 
'^ieyiz ^'"s correlated in a very complex manner. 



and 0^f = (p^'^'iEl 



nz(ie) 
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It should be 



^1 The iz = 1 bin corresponds to the transit of the source. 

^' Each run lasts ~30niin, OFF-source runs being used to estimate the hadronic background. 
'^ The average is taken over the bin according to the physical meaning of the flux, i.e. by 
integrating over E or log(i<^) in a dN/dE or in a i^F^(= E"^ x dN/dE) representation, respectively. 



CAT OBSERVATIONS OF MRK 501 AND MRK 421 
Spectral properties of Mrk 501 in 1997 and 1998 
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FIGURE 1. (a) Mrk 501 nightly integral flux above 250 GeV in 1997 and 1998; (b) Mrk 501 
VHE SED between 330 GeV and 13 TeV for three independent data-subsets, corresponding to 
different activity states in 1997 (see [4] for details). 



Mrk 501 exhibited a remarkable series of flares during the whole year 1997, going 
down to a much lower mean flux in 1998 (Fig. la). The detailed spectral analysis 
of the 1997 data can be found in [4]; here we extend it to the 1998 results. 

The VHE spectral energy distribution (SED) of Mrk 501, derived for different 
flaring- activity states in 1997 (Fig. lb), shows a significant curvature which is now 
well confirmed by different ground-based experiments [1,6]. The peak 7-energy is 
found to lie just above the CAT threshold, and it seems to shift towards higher 
energies as the flux increases. To check this spectral variability by a more robust 
method, the hardness ratio has been computed for five different-level intensities: 
the correlation observed in Fig. 2a confirms the hardening of the VHE SED during 
flaring periods. Fig. 2b shows the broad-band SED of Mrk 501 for two flaring 
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FIGURE 2. (a) Hardness-ratio (HR - „ 

^ ^ ^ -'Ve>450 GoV ' 

10"" cni~2g-i). ^5j ]^j.]j 501 X-ray and VHE spectra for April 7"> and 16**^. The EGRET upper 
limit corresponds to observations between April 9**^ and 15**^. 



dates in April 1997: it exhibits the two-bump structure typical of blazars, with a 
dip indicated by the contemporary EGRET upper-limit point in the GeV energy 
range [11]. The obvious correlation of the X-ray emission from BeppoSAX data [8] 
with TeV emission strongly suggests that the same particle population is responsible 
for emission in both energy-ranges, i.e. it supports the picture given by leptonic 
models [5], in which an energetic electron beam propagating in the magnetized 
plasma jet produces X-rays through synchrotron radiation as well as VHE 7-rays 
through inverse Compton scattering of low-energy photons. 
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FIGURE 3. (a) Mrk 501 VHE SED from 330 GeV to 5.2 TeV in 1998, as compared to the MF 
data subset in 1997 (see Table 1); (b) 68%, 90% and 95% CL contours in the {7'=^,/3'=^} plane 
for the 1997 AV and 1998 data sets (see Table 1); each contour is obtained by projecting the 
3-dimension ellipsoid along the (J)q^ axis. 

Unlike the 1997 spectra, the spectrum of Mrk 501 in 1998 does not show any 
curvature (Fig. 3a), indicating a VHE peak energy well below the CAT threshold at 
that time. Fig. 3b confirms in particular that, in spite of larger statistical errors, the 
1998 data cannot be fitted by the average spectrum shape found in 1997 (AV data 
set, see Table 1), as would be expected in the absence of any spectral variability. 
With a mean flux much lower than that of 1997, the power-law shape found in 1998 
is therefore consistent with the previously suggested scenario, in which the peak 
7-energy seems to be correlated with the VHE flux. 



Mrk 421's spectrum in 1998 



As reported in [10], Mrk 421 is the second extragalactic source detected by CAT: 
almost quiet in 1996-97, the source showed small bursts in 1997-98, together with 
a higher mean activity (Fig. 4a). The energy spectrum derived for the 1998 flaring 
periods is well represented by a simple power law (Fig. 4b) with a differential 
spectral index 7 = 2.96 ± 0.13***'^* ± 0.05^^*^* [10]. This agrees with a recent result 
of the HEGRA group concerning the 1997-98 period [2]. It is however in contrast 
with the former behaviour of the source at the time of the 1995-96 flaring period, as 
observed by the Whipple group, who found 7 = 2.54 ± 0.03^**^* ± 0.10*^^*^ [6]. Thus, 



; Dec 


Jon 


Feb 


Mor 


Apr 


Uo, 


L 1996-97 












r 




■■ 










,1 




t||. 


,1 t 


-f-- 


'""1 


,' 


""i 




"""Trt""" 


L fl 


1 1 


1,, 


\~} 


IT 





50420 50440 50460 50480 50500 50520 50540 50560 50580 



> 

H 



; Dec Jon 
L 1997-98 1 

1 


-if- 


Feb 


Mar 


1 
-<^#- 


Apr 


May 








t 




T T 









50800 50S20 50840 50S60 50880 50900 50920 50940 

Modified Julian Day 




Energy (TeV) 



FIGURE 4. (a) Mrk 421 nightly integral flux above 250 GeV between December 1996 and May 
1998. The dashed line represents the mean flux over the two years; (b) Differential flux of Mrk 421 
for the flaring periods in 1998. Only the energy bins shown with fuUfilled circles (from 330 GeV 
to 5.2 TeV) were used in the likelihood method (see text). 

the higher value of the differential spectral index found in this work could come 
from a correlation between the intensity level and the spectral hardness, like that 
observed for Mrk 501 in 1997 (see [4] and above). A similar evidence for Mrk 421 
was in fact suggested in [13] on the basis of the Whipple 1995-96 low-flux data. 

In the framework of leptonic models [5], which succesfully explain the Mrk 501 
broad-band SED in 1997, the absence of any obvious spectral curvature reported 
by all experiments implies in any case that the peak energy of the inverse-Compton 
contribution to Mrk 421's spectrum is significantly lower than the CAT detection 
threshold, as for Mrk 501 in 1998. This is not surprising since the corresponding 
synchrotron peak is lower than that of Mrk 501 in 1997, and since we have seen 
that leptonic models predict a strong correlation between X-rays and 7-rays. In 
fact, such a correlation was directly observed on Mrk 421 in Spring 1998, during 
a coordinated observation campaign involving ground-based Cherenkov imaging 
telescopes (Whipple, HEGRA, and CAT) and the ASCA X-ray sateUite [12]. 



CONCLUSION 



Since Mrk 501 and Mrk 421 lie at the the same redshift (~0.03), spectral differ- 
ences between them must be intrinsic and not due, in particular, to absorption by 
the diffuse infrared background radiation. This allows direct and relevant compari- 
son of their spectral properties. The correlation between X-ray and 7-ray emissions 
is now proven for both sources, supporting the simple and most natural scenario 
given by leptonic models [5] for the origin of blazar TeV flares, in which a single 
leptonic population is injected into the radio jets and produces correlated X-ray 
synchrotron and VHE 7-ray inverse Compton radiations. CAT observations com- 
plete this picture with some evidence of spectral variability of Mrk 501 in 1997, 
suggesting that this leptonic population is responsible for the hardening of the en- 



TABLE 1. Best fitted spectral parameters of Mrk 501 and Mrk 421 in 1997 and 1998 
from the T^p' and 'W^ shape assumptions. The data-subsets of Mrk 501 in 1997 are 
defined in Fig. lb, while the average set (AV) contains the data of the whole year. Fluxes 
are given in units of 10"^"'^ cm~^s~^TeV~^. If, according to A, the T-f^^ hypothesis is 
favored, then the last two columns gives the curvature term [3 and the peak-emission 
energy E^^ = 10^^. If not, /3 is quoted between brackets. 
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ppcak 
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Mrk 501'97 LF 
MF 
HF 

AV 


13.6 

40.5 

3.1 

57.2 


10.7 
47.1 
29.1 
61.5 


3.13 ±0.19 

4.72 ±0.14 

17.60 ±0.61 

5.19 ±0.13 


2.32 ±0.09 
2.25 ±0.05 
2.07 ±0.04 
2.24 ±0.04 


0.41 ±0.17 
0.52 ±0.08 
0.45 ±0.09 
0.50 ±0.07 


410 ±201 
583 ± 104 
840 ± 108 

578 ± 98 


Mrk 501 '98 


3.3 


0.09 


1.25 ±0.16 


2.97 ±0.20 


(0.21 ±0.73) 


- 


Mrk 421 '98 


5.1 


0.34 


1.96 ±0.20 


2.96 ±0.13 


(0.28 ±0.49) 


- 



^ total time (h) ON-source. 

tire high-energy part of the electromagnetic spectrum during flares. To date, the 
peak 7-energy of Mrk 421 has always remained well below the CAT threshold, pre- 
cluding a more accurate spectral study. Therefore, testing VHE spectral variability 
as a general feature of blazars requires more multiwavelength observations with 
a large dynamic range in intensity. In any case, the comparison of the spectral 
properties of Mrk 501 and Mrk 421 in 1997 and 1998 fits the picture in which the 
VHE peak 7-energy shifts with increasing VHE fiux: as shown in Table 1, spectral 
curvature seems to be characteristic of high fiaring-activity states (i.e. Mrk 501 in 
1997), whereas low-activity spectra (i.e. Mrk 501 and Mrk 421 in 1998) are always 
compatible with pure power-laws. 
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